Parameterized Beyond-Einstein Growth 
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A single parameter, the gravitational growth index 7, succeeds in characterizing the growth of 
density perturbations in the linear regime separately from the effects of the cosmic expansion. The 
parameter is restricted to a very narrow range for models of dark energy obeying the laws of general 
relativity but can take on distinctly different values in models of beyond- Einstein gravity. Motivated 
by the parameterized post-Newtonian (PPN) formalism for testing gravity, we analytically derive 
and extend the gravitational growth index, or Minimal Modified Gravity, approach to parameterizing 
beyond-Einstein cosmology. The analytic formalism demonstrates how to apply the growth index 
parameter to early dark energy, time-varying gravity, DGP braneworld gravity, and some scalar- 
tensor gravity. 



I. INTRODUCTION 

The acceleration of the cosmic expansion points to new 
physics beyond the standard models of particle physics 
or gravitation. But the nature of this physics is not 
clear. While current data constraints are consistent with 
Einstein's cosmological constant A, the uncertainties are 
still substantial. Even the simplest deviation from this 
picture, an assumed constant effective equation of state 
w, is determined to no better than 20% (at 95% confi- 
dence level) when combining the full set of current super- 
nova distance, cosmic microwave background, and baryon 
acoustic oscillation measurements, with systematic un- 
certainties included Given that virtually no physical 
explanation beyond A predicts a simple, purely constant 
equation of state, it is clear that we cannot claim to have 
zeroed in on Einsteinian physics as the solution. 

As we can look beyond Einstein's cosmological con- 
stant, so too we can look beyond the framework of Ein- 
stein gravity. To address the question of the nature of the 
new physics requires measuring both the expansion his- 
tory of the universe and the history of structure growth. 
The growth history can provide independent information 
on the accelerating physics in the case when gravity de- 
viates from general relativity (e.g. [1, 0, 0, @|). As we 
describe in detail below, the structure growth history is 
conveniently described by a parameter 7, first introduced 
by [|| and expanded upon by 0] , within Einstein gravity, 
to supplement the conventional cosmological parameters 
of the matter density and the equation of state of the 
dark energy. 

The goal of this article is to motivate and explore fur- 
ther the idea of a model independent parametrization of 
the deviations of growth from Einstein gravity. Just as 
the model independent parametrization of the expansion 
history equation of state w(z) proves valuable even with- 
out a particular Lagrangian in mind, we use 7 without the 
intent to adopt a specific gravitational action, though we 
do show results for some specific (toy) models. We find 
that 7 is confined to a narrow range of values near 6/11, 
independent of the nature of the dark energy if general 
relativity holds. In contrast, beyond-Einstein models of 
gravity can lead to rather different values for 7. This 



difference then establishes a standard for measurements 
of the growth of structure required to exclude alternate 
models of gravity. 

Linder Q proposed that this independent information, 
apart from the expansion physics, be parameterized by a 
single gravitational growth index 7, defined through 



, , f"(da'/a') [O m (a') T -l] 

/(a) = e J o J , 



(1) 



where J7 m (a) is the matter density fraction of the total 
energy density as a function of scale factor, g — S/a and 
5 = 5pm/ Pm is the linear density perturbation in the mat- 
ter. (See 0] for an early approach to the growth index 
within general relativity and tied to the expansion his- 
tory.) In [9j, this was combined with the expansion his- 
tory into a full cosmological fitting framework called Min- 
imal Modified Gravity (MMG), using that the expansion 
history of the universe, a(t), can be phrased in terms of 
an effective equation of state ratio w(a) = wq + w a (1 — a) 
shown to be extremely successful for describing a wide 
variety of physics, including gravitational modifications 
[l(J[llj]. Comparison of the parameters {wq, w a , 7} fitted 
to data (see examples in [9J) could establish statistically 
significant evidence for beyond-Einstein physics and pro- 
vide information on its origin. 

In <JTT]we reexamine and derive analytically certain im- 
portant properties of the parameter 7, which were origi- 
nally obtained numerically. The physical motivation for 
the growth index approximation to the growth history 
is strengthened in mill and evaluated for time varying 
and early dark energy. We derive a relationship between 
the growth index and the PPN formalism for gravity 
theories in mVl and show the success of the growth in- 
dex in characterizing DGP braneworld gravity and some 
scalar-tensor gravity theories. The analytic derivations 
serve an important role as a foundation for understanding 
both extensions and breakdowns of the growth index. Of 
course exact calculations of growth within specific gravi- 
tational models, not treated here, play an important role 
as well. 
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II. PARAMETERIZING GROWTH 

Considering linear density perturbations in the matter, 
8 = &Pml Prm the equation for their growth within general 
relativity is given by 



d 2 S 



dS 



— + 2H(a)- -4^ 



0, 



(2) 



where H = a/ a is the Hubble parameter. This can be 
rewritten as 



dG 
din a 



ldlniP 
2 dhia 
ldlnH 2 
2 d In a 



G + G 



0, 



(3) 



where G = d\n{8 / a) / d\na. Note that by normalizing 5 
by a we remove the pure matter universe growth behav- 
ior, which would give G — 0. The dimensionless matter 
density is S! m (a) = Q m a~ 3 /[H(a)/H ] 2 , where f2 m is its 
present value. 

The growth history can be solved formally by quadra- 
ture, giving the solution 



G(a) = -1 + [a^Hia)}- 1 /° ^ a'^Hia') 

Jo a' 



l + ^ m (a')-G 2 («') 



(4) 



This can also be written in terms of the dark energy 
density fl w (a) = 1 — f2 m (a), for a flat universe, as 



G(a) = -1 + 



Jo a 
b --' d -n w {a')-G 2 {a') 



(5) 



For growth during the matter-dominated era, G will be 
small, and this holds reasonably well even as dark energy 
comes to dominate, since fl m is still non-negligible today: 
for the concordance cosmology, today G is of order —1/2. 
So let us neglect the G 2 term in Eqs. (HI)-© - while this 
is motivated by the asymptotic growth behavior during 
matter domination, it should be a reasonable approxima- 
tion throughout the growth history. In practice we find 
that g(a) — S/a agrees with the exact solution to ~ 0.2%. 
This now gives an explicit solution for the growth. 

We can use 



H 2 /H 2 = n m a- 3 [l + n w {a)/n m (a)} 



(6) 



to write at the same level of approximation (that the 
universe is not too far from matter-dominated), 



G(a) 



1 



-Cl w (a) 



-5/2 



a ^ a <^n w (a'). (7) 
o a 



For any particular model of H(a) 7 or Q. m (a) or £l w (a), 
we can then evaluate the growth history. 



Connecting this new analytic approach to the growth 
index defined in Eq. {!]), we find immediately 



G(a) = n m ( a )-y - 1. 



(8) 



In the same approximation used to derive the quadra- 
tures for G, we can write 



jK-G/n w {a), 
and then using Eq. (|7|) we have 

2 4 Jo a 

1 1 f 1 du 



(9) 



2 4 ,/ u 



u 5 / 2 fl w (au)/Q w (a). 



(10) 



To first order in deviations from matter domination, 
e.g. at early times, fl w (a) ~ a~ 3u, °°, where = w(a 
1), and we can evaluate Eq. (fTTJ|) to obtain 



7o 



3(1 - Wqq) 

5 - 6^00 



(11) 



Equivalently, we can substitute Eq. into Eq. ([3]) and 
find 



7o 



3 + G 
5 + 2G" 



(12) 



where G = dhiG / d\na = — 3woo to this level of approx- 
imation. 

These asymptotic values accord exceedingly well with 
the fitting formula of 0] . We can write 



7o 



6 -3(1 + woo) _ 6 3 

+ T7Tt(1 + w oo)- (13) 



11 -6(1 + k^) 11 121 



We can compare this to the fitting formula from [8( that 
gives 

7 = 0.55 + 0.02 [l + w(z = 1)], w<-l (14) 
= 0.55 + 0.05 [l + w(z = 1)], w > -1 (15) 

over the whole range < a < 1. For the cosmological 
constant case, w = —1, the asymptotic value is joo = 
6/11 = 0.545 compared to the numerically obtained fit 
7 = 0.55 for the whole growth history. 

The first order correction to the cosmological constant 
case for a different value of constant w is quite small 
for the asymptotic formula (fl~3|) . with a coefficient of 
0.025 times 1 + w. This agrees with the fitting correc- 
tion from Eqs. (p~4|) - (p~5|) . which have a coefficient of 0.02 
(0.05) times 1 + w for the case when 1 + w < (>) 0. 
The asymptotic value is closer to the fit for the w < — 1 
case because then matter domination lasts longer. The 
asymptotic form is also in agreement with the pioneering 
work of 0|, who expanded the growth equation about 
^m(a) = 1 for constant w. Their asymptotic term agrees 
with Eq. (jlip and they find the next order is given by 
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(15/1331) [1 — ft m (a)) for the w = — 1 case; even when 
1 — f2 m (a) is much larger than the ~ 10~ 9 for ACDM at 
CMB last scattering, this correction is negligible. 

The new elements in Eq. (fTU|) include the integral na- 
ture of the relation, allowing for treatment of the whole 
growth history to some redshift, rather than an instan- 
taneous measure of growth (thus allowing a prediction 
of Gtoday unlike previous work), the ability to treat time 
varying dark energy, and the clear identification of the 
key assumptions such as early matter domination that 
allow the single parameter to succeed. 

We emphasize that the numerical fitting form of [8| 
covers the entire growth history, not just the asymptotic 
high redshift behavior, and does not assume Gel but 
rather fits to the exact numerical solution. Nevertheless, 
it is instructive to pursue the analytic arguments further 
to understand the motivation for phrasing the growth 
history in terms of a gravitational growth index, and why 
a single parameter proves so successful, even in the case 
of time varying equations of state. 



III. GROWTH AND DARK ENERGY 

One question to ask is why the fitting form ([1]) is phys- 
ically appropriate. In the asymptotic limit, G will be lin- 
early proportional to il w (a) = 1 — f2 TO (a) since Eq. ([3]) 
is linear (for small G; cf. Eq. [7|), but one can imagine 
other forms besides fi m (a) 7 — 1 that have this property. 
We consider three possibilities for fitting forms of G: the 
standard one of f2 m (a) 7 — 1, one directly proportional to 
fl w (a), and lnf2 m (a) 7 . Each of these has the appropri- 
ate limit that G vanishes linearly in 1 — f2 m (a) as a — > 0. 
Defining 7 for the entire growth history in terms of these 
forms gives the options 



7 



_ ln(l + G) 



G 



G 



]nil m (a) ' Q w (a) lnO m (a) 



(16) 



As a guide to defining 7 as a useful parameter (e.g. 
nearly constant) over the entire growth history from 
a = to the present, we can examine its late time be- 
havior, when o» 1. For dark energy domination in the 
future, the matter density f£ m (a) — > and matter density 
perturbations S stop growing, so G — » — 1, in agreement 
with Eq. (@|. The second definition of 7 in Eq. (fTB)) 
leads to 7 — » 1 in the future; the third definition has 
7^0. Only the first, original definition of the growth 
index preserves its stability over the entire growth his- 
tory from asymptotic past to asymptotic future. Indeed, 
7 varies from its present value by less than 2% back to 
a -C 1 and 6% forward to a = 5, for the cosmological 
constant case. Since 7 enters into an integral relation, 
then even this mild variation is further smoothed over, 
giving a constant 7 as an excellent approximation. Q 
found it to reproduce the exact growth for a wide variety 
of models to better than 0.2%. Overall, the definition of 
a gravitational growth index through Eq. ((T|) is therefore 
physically well motivated. 



We can carry this further, showing that a single growth 
parameter suffices rather than needing a function (of red- 
shift). Considering dynamical dark energy, rather than 
a constant w model, expand 



Q w (a 



(17) 



in Eq. (|10[) . where we approximate the time variation 
of the equation of state at high redshift as a power law 
correction to the asymptotic Woq value. The solution for 
the growth index is 



7 



3(1 - Wqq) 

5 - 6ui„o 



Bx a x 



(5 - 6woo)(5 - 6woo + 2x) 



(18) 



For the standard equation of state parameterization 
w(a) = wq + w a (l — a) = Woo — w a a, we have B = 3w a , 
x = 1 . The correction to 7 is exceedingly small, at early 
times proportional to a, with furthermore a small co- 
efficient -(3/143)w a , generally < 0.01. So even in the 
dynamical dark energy case, 7 is found to be nearly con- 
stant. (The growth index fit was examined numerically 
in @.) 

Suppose we consider early dark energy, where its en- 
ergy density is not completely negligible at high redshifts, 
but possibly contributing up to a couple percent of the 
total energy density at CMB last scattering [12, 0, [H| • 
We might expect this to upset the gravitational growth 
index since early dark energy directly affects the matter 
domination used to derive the analytic asymptotic behav- 
ior. In order for dark energy to contribute non-negligibly 
to the early energy density, its equation of state must ap- 
proach w = 0. Substituting w M = in Eq. (fTTj) . we see 
7oo = 3/5. However, we can actually obtain this solution 
without approximating G < 1; the exact solution for 
growth in a universe with two components, each evolv- 
ing (at early times) with density proportional to a -3 , as 
ordinar y m atter, but with a fraction Q e not clustering, is 



ordinar y m 



9 = 5/a 



(-5+V25-24n e )/4 ^ -(3/5)0, 



(19) 



or 7 = 3/5. (To next order, 7 « 3/5 + (3/125)Q e , in 
agreement with Q.) 

The redshift distortion factor (3 [l?} used in galaxy red- 
shift surveys can be written as (3 ~ 1 + G = f^, and one 
often sees il^f used. For the concordance model we have 
seen that 7 = 0.55 is a much better fit. Using il^f intro- 
duces a needless 6% error in (3, or a systematic bias of 
0.03 (11%) in the value of fi m derived, so more accurate 
results will be achieved if f^ 55 is used in place of fi„ 6 . 



IV. GROWTH AND GRAVITY 

To this point, the growth behavior has been completely 
specified by the expansion history. That is, in Eq. ([3]) 
the terms only involve H(a) or O m (a), and so the grav- 
itational growth index 7 can be defined in terms of the 
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effective equation of state w(a). Furthermore, we have 
seen that 7 has quite a small range over a reasonably 
large variety of dark energy models, i.e. for w = —1 to 
— 1/3 (including effective curvature energy), 7 only varies 
over 0.55-0.57. 

Now we consider what happens when we allow alter- 
ations to the gravity theory. There is no unique prescrip- 
tion for how modified gravity theories affect the growth 
equation, even in the linear regime, though see [8( for 
an attempt to provide a somewhat general, if formal, 
treatment. Effects on the growth equation include the 
introduction of scale dependence, anisotropic stress, and 
variation of the gravitational coupling (i.e. Newton's con- 
stant). We briefly discuss the first two of these in ^IV Dl 
and [V] and concentrate here on the last one. Again 
we emphasize that our goal is a model independent ap- 
proach, rather than adopting a specific theory. 



A. Varying gravity 



effective energy density il w (a). We discuss this further 
in the specific DGP braneworld and scalar-tensor gravity 
examples below. 

In the scaling case, we can write Q — 1 = Afl w (a) and 
find the asymptotic growth index to be 



7oo 



3(1 



A) 



5 — 6w c 



(23) 



We have verified that, as in ^III[ it is best to define 7 
in the usual way through G = f2 m (a) 7 — 1, not G — 
[Qf2 m (a)] 7 — 1, even in the presence of beyond-Einstein 
gravity leading to Q / 1. That is, the gravitational de- 
viation Q — 1 enters strictly through the growth index 7, 
while the expansion history determines the growth his- 
tory G apart from the value of 7. This separation of 
physical effects from the expansion rate and from the 
gravity theory through distinct beyond-Einstein param- 
eterizations w and 7 is an important point that clarifies 
the nature of the beyond-Einstein physics. 



The last term of Eq. © contains a hidden dependence 
on the gravitational coupling, multiplying the source 
term Q m (a). Here we rewrite the growth equation, ex- 
plicitly showing the effect: 



dG 
din a 



, 1 d:nir- \ ri . 
4+ 2^hTT )G + G 



1 dlnH 2 

2 am a 



-[l + (Q(o)-l)]n m (a)=0, 



(20) 



where Q — 1 gives the fractional deviation of the coupling 
from the general relativistic case (i.e. Q — 1 means the 
coupling is given by Newton's constant). 

We can now repeat the analysis of £|TT1 finding 



G{a) = -l+(a i H)- 1 C — a'^H 
Jo a 1 



1 + -Q(a')fl m (a') 



This is turn leads to a revised gravitational growth index 



7* o 



' — u 5 / 2 fl w (au)/Q w (a) 



1 du 



u 5/2 [Q{au) - l]/n w (a). (22) 



There are three cases to consider for the gravitational 
deviation. If (Q — 1) ~ a q at early times, and q > Sw^, 
then the modification to 7 from the altered gravitational 
coupling will be negligible. However, if q < — 3u>oo then 
there will be strong modification. In fact, in this case the 
usual matter-dominated growth behavior at high redshift 
is broken (G = 0, or S ~ a, is no longer a solution at early 
times), and large scale structure would not accord with 
observations. This leaves the main, and physically best 
motivated, case of the scaling behavior q = Sw^ where 
the same physics responsible for the variation in the force 
of gravity also affects the expansion rate, giving rise to an 



B. Braneworld gravity 

To test our new expression for the gravitational growth 
index we consider several examples of gravity beyond 
general relativity. First, we examine a theory altering 
the Einstein-Hilbert action by a term arising from large 
extra dimensions, the DGP braneworld theory (TR 
On scales smaller than the Hubble scale but still within 
the linear density perturbation regime, the effect on the 
growth behavior is that of a variation in gravitational 
coupling, with [2| 



(q-i; 



DGP 



1 A-^(a) 



3 \l + Cll(a) 



(24) 



(21) 1) 



Note that, as predicted, one has the scaling behavior (Q- 



asymptotically and 



.4 



Q-i _ 1 i + n m (a) 

l-n m (o)~ 3 1 + f! 2 m (a) 



(25) 



is of order unity. This ensures that the gravitational de- 
viation neither violates matter domination nor has neg- 
ligible effect on growth. 

The effective equation of state from DGP gravity is 
w(a) — — 1/[1 + fl m (a)] 0], approaching Woo = —1/2 at 
high redshift. Substituting the DGP values for A and 
Wqo m to Eq. (|23[) . we find the growth index 



7oo,DGP 



11 
16 



0.6875. 



(26) 



This accords exactly with the asymptotic numerical so- 
lution and extremely well with the gravitational growth 
index fit 7 = 0.68 over the whole history, given by 0]. 

If for any point in the growth history we naively sub- 
stitute into Eq. (|2"3"|) the values of A from Eq. (f2"S"|) and 
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w(a) = — 1/[1 + f2 m (a)], we obtain 

„ 7 + 5r> m ( a ) + 7», 2 T1 ( a ) + 3^( a ) 
7DGP ~ [l + ^(a)][ll + 5fi m (a)] ' 1 ' 

which agrees well with the exact numerical solution. The 
growth index takes on the value 11/16 = 0.6875 in the 
asymptotic past and 0.634 in the asymptotic future (vs. 
7/11 = 0.636 from the formula, which was derived under 
matter domination). At the present the numerical value 
is 0.665 (for Cl m = 0.28, vs. 0.674 from the formula), so 
the single parameter 7 = 0.68 is an excellent approxima- 
tion for growth at any time through the present, holding 
constant to 2%. In addition to the remarkable constancy 
of the growth index, 7 stays well distinct of the pure 
expansion history prediction within general relativity of 
7 = 0.55 — 0.56 as the DGP equation of state evolves 
from w = — 1/2 in the past to —1 in the future. 

We illustrate this important property of separation of 
growth history parametrization from expansion history 
in Figure [T] Here the growth index is shown as a ratio to 
the (exact) general relativity (GR) value for the cases of 
quintessence and braneworld models. For the braneworld 
case, being a single parameter model, w(z = 1) deter- 
mines Q m (note f2 m will be far from 0.28 when w{z = 1) 
is far from —0.62). We see that the curve of 7 as we vary 
the expansion history parameter w{z = 1) is quite flat, 
showing success in obtaining a growth parameter distinct 
from expansion effects. 

For the quintessence case (taking Q m = 0.28, u> = — 1, 
and then w(z — 1) serves as a proxy for w a ) the curve 
representing the fitting form Eq. ([T5]) (with a coefficient 
0.04 rather than 0.05 as a better fit over the restricted 
range w(z — 1) <E [—1,-0.5]) is within 0.2% of unity, 
showing the success of this fitting form. 

The growth index formalism thus possesses these im- 
portant properties: 1) the constancy of 7 in redshift, 
allowing a single parameter description of the gravity de- 
viations beyond-Einstein, 2) the independence of 7 from 
the expansion history, i.e. separating out the expansion 
influence on the growth so as to give a distinct window 
on the gravitational physics, and 3) clear signal of the 
origin of the beyond-Einstein physics, achieving in the 
braneworld case over 20% deviation from the general rel- 
ativity prediction - while the "noise" from the expansion 
history w within general relativity affects 7 at the 0.2% 
level, a factor of 100 less. Thus, the gravitational growth 
index provides a clear and effective parameterization of 
beyond-Einstein gravity. (While the theory signal-to- 
noise is high, achieving observational constraints is more 
challenging, with [9( estimating that next generation ex- 
periments will determine 7 to within 8%.) 



C. Scalar-tensor gravity 

Next we consider scalar-tensor theories of gravity, in- 
volving coupling of a scalar field 4> to the Ricci curva- 
ture R of the form of F(<j>) R in the action. Because of 
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FIG. 1: In the parameterized beyond-Einstein approach, the 
growth index 7 tests gravity separately from the expansion 
and clearly distinguishes beyond-Einstein gravity from physi- 
cal dark energy. The curves show the values of 7 change little 
as the expansion history, here represented by the equation of 
state w(z = 1), varies. Solid curves show the ratio of the nu- 
merical solution for 7 in the braneworld case to that in general 
relativity (GR), and in the quintessence case the ratio of the 
fitting form 7 = 0.55 + 0.04 [1 + w(z — 1)] to the numerical 
GR solution. (Dotted curves assume 7 = 0.55, the GR value 
in a cosmological constant model.) 



the physical coupling between the modification of the ex- 
pansion and the modification of the growth equation, we 
again might expect the scaling relation to hold, where the 
gravitational deviation Q — 1 is of order the effective en- 
ergy density Q, w (a) in the expansion, at least for theories 
consistent with observations of large scale structure. This 
would give an appreciable, but not pathological, influence 
on growth. Scalar-tensor theories where the scalar field 
is responsible for the current acceleration are often called 
extended quintessence theories [20| . 

In extended quintessence (EQ), one can show that the 
main effects on expansion history come from the scalar 
potential (self-interaction) and the change in the gravita- 
tional coupling (SttGn)^ 1 —> F, where Gn is Newton's 
constant. For the growth equation, deviations arise pri- 
marily from the variation of the coupling, with additional 
terms being suppressed by factors of (k/H)~ 2 , where k 
is the wavemode of the density perturbation. For calcu- 
lations of w(a) and the modified Poisson equation, see 

m. 

The effective modification to the growth history can 
be treated through a deviation in the source term, with 
Q — 1 = (SttGnF)^ 1 — 1. Given a form for F, this can be 
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used in Eq. (|23|) to obtain the gravitational growth index 
for the scalar-tensor theory. In EQ, the i?-boost mecha- 
nism [HJ operating during matter domination drives the 
theory toward general relativity, so we have a consistent 
picture of the usual matter-dominated growth around the 
time of CMB last scattering. As the scalar field comes to 
dominate in the late universe, the theory may diverge 
from Einstein gravity. (Because we are exploring the 
growth index formalism, we do not here worry about 
constraints from solar system tests of gravity or about 
specific scalar-tensor theories.) 

To give some flavor of calculating 7, we adopt a toy 
model with coupling 



F 



SttGn 



1 



B 



1 + (a/a*)-« 



(28) 



where B is the amplitude of the coupling variation and 
a* is a transition scale factor. At early times, (Q — 1) ~ 
(a/a*) 9 . The dark energy density (the quintessence part 



of EQ) varies as fl w (a) 



-3w 



00 . As before, if q > 



—SiUoo, the growth deviation arising from (Q — l)/O w (a) 
in Eq. (|23p will be negligible; if q < — 3u>oo then the 
growth source term will be drastically affected at early 
times. As mentioned before, the coupling between the 
scalar field evolution, and hence gravity deviation, and 
the expansion provides a motivation for the scaling be- 
havior q = —3 Wqq. In this regime, 



.4 



1 



Q w (a) 



B 



,3iUo 



1 - Or, 



(29) 



Substituting this into Eq. (f2"3"| predicts 7 = 0.571 for 
a* = 0.5, B = 0.03, and Woo = — 1. Numerical solution 
of the growth equation gives 7 = 0.571 at high redshift, 
7 = 0.564 today. This is close to the Einstein range since 
B is small, but the main point is that 7 is quite constant 
over the growth history. However, for parameter val- 
ues such that the gravity deviation causes A to approach 
unity, 7 does start to vary (but of course there would 
be severe departures from general relativity today). As 
mentioned, this was a toy model giving but a brief taste 
of the rich phenomenology of scalar-tensor gravity and 
a full analysis should take into account all the observa- 
tional constraints. (Also see [23| for analysis of linear 
perturbations within f(R) gravity.) 



D. Relation to PPN 

As we survey a larger range of gravitational theo- 
ries, the situation becomes more complex. The matter 
source term in the growth equation does not arise purely 
from the Newton-Poisson equation V 2 $at = AirSp, or 
<I>Ar = —4irk~ 2 pa 3 (5/a) in Fourier space, where <i>jv is 
the Newtonian potential. The equations of motion actu- 
ally depend on two potentials, $ and appearing in the 
metric as 



(written for a flat universe and in longitudinal gauge for 
simplicity) . 

The second order differential equation for the density 
contrast, a la Eq. @, comes from both the density and 
velocity perturbation evolution, and involves both ^ and 
$. The combination 



* + $ = -fc- 2 n, 



(31) 



is not necessarily zero, as in general relativity, in the 
presence of anisotropic stress II. While saying anything 
general about growth is difficult, let us motivate one ap- 
proach. 

The fractional correction to the source term was de- 
fined as Q — 1 in £|IV Al For wavemodes on sub-Hubble 
scales, the dominant correction to the growth equation 
takes the form (following the derivations of 0, [Hj]) 



k 2 ^ 



-fc 2 * + (2/3) fc 4 (* + $)/p„ 



(32) 



ds 2 



-(1 + 2^)dt 2 + a 2 (t) (1 + 2$)df 2 , 



(30) 



We note that we expect the scaling behavior to hold be- 
cause the anisotropic stress is generated by the dark en- 
ergy fluid. The deviation Q — 1, and the factor A in 
Eq. (f2"3"|) , is proportional to l+Q/ty. We can now consider 
a direct connection to the parameterized post-Newtonian 
(PPN) formalism (26j. The first PPN parameter is just 
the ratio of the first order corrections of the gu and goo 
parts of the metric, hence equivalent to — $/^. This 
PPN parameter is unfortunately also given the symbol 
7, so we will denote it as 7ppn- Following £ )IV Al we see 
a correction to the growth index, 

= 3(l- Woo + (2/9)(fc/H) 2 [l- 7P PN]) (33) 
5 - Qwoo 

Note that in this case 7 has a scale dependence. (For 
large wavemodes k, or small wavelengths, the k 2 behav- 
ior is cut off by other terms we neglected, such as the 
usual Jeans term of the matter pressure.) As the beyond- 
Einstein deviation 1 — 7ppn vanishes, the growth index 
approaches the general relativity expression. This illus- 
trates but the barest outline of the complexity of beyond- 
Einstein gravity. 



V. CONCLUSION 

Observations, while broadly consistent with the con- 
cordance cosmology of Einstein's cosmological constant 
within Einstein's general relativity, still offer significant 
leeway for beyond-Einstein physics. Whether there is 
such a deviation and how to distinguish its physical ori- 
gin - from a new field or an extended theory of gravity 
- are major questions to address with the next genera- 
tion of experiments. However on the theory side, there 
is no universally accepted model beyond-Einstein. This 
drives us to develop a usable, model independent formal- 
ism, along the lines of the parameterized post-Newtonian 
approach, in which to evaluate future data and guide us 
to a theory. 
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Here we have demonstrated that the Minimal Modi- 
fied Gravity (MMG) approach of combining a gravita- 
tional growth index 7 to describe the mass perturbation 
growth history with equation of state w(a) parameteri- 
zation to describe the expansion history is in many cases 
robust, accurate, and broadly applicable. The growth in- 
dex 7 clarifies the nature of beyond-Einstein physics by 
separating the expansion effects on the growth from the 
gravitational theory effects. 

Rather than needing a full function, the single growth 
index parameter 7 (a la PPN) is extraordinarily success- 
ful: 1) it reproduces the exact growth history in a highly 
accurate manner (< 0.2% deviation), 2) stays constant 
to high accuracy, 3) describes both physical dark en- 
ergy (iTTTJ) and beyond-Einstein gravity models consis- 
tent with observations f ^IVI) . and 4) exhibits clear dis- 
tinction between these different physical origins by de- 
viating as the gravity theory does (up to ~ 20% from 
general relativity) even when the expansion histories are 
identical. 

We emphasize, however, that while MMG is surpris- 
ingly broad in its physical realism and parameterized 
beyond-Einstein growth is an intriguing formalism for 
exploring new physics, there is considerable work still 
to do. This article has dealt only with the linear den- 
sity perturbation regime of growth for modes below the 
Hubble scale; other cases, in particular the nonlinear 
regime, may not be susceptible to a model independent 
approach. However, [9| suggests first steps for addressing 
the translinear regime, 6 ~ 1, within MMG; this region 



will be important for weak gravitational lensing probes 
of cosmology. Modes approaching the Hubble scale can 
be probed to some extent by the integrated Sachs- Wolfe 
effect, and these observations may offer some clue to 
beyond-Einstein physics (see, e.g., [27]). The Hubble 
scale, and other features like a Vainshtein or Yukawa 
scale, can introduce scale dependence in the growth as 
well [281 [29( . Microphysics such as a dark energy sound 
speed or anisotropic stress (or coupling matter to dark 
energy, or non- minimal coupling to gravity), can confuse 
the interpretation of the growth history, as outlined by 
3C I . A specific formal model for this was proposed by 
3l| . though the dark energy perturbations become large, 
possibly giving observational difficulties. 

Our understanding of the full range of models beyond- 
Einstein has a long way to go, but Minimal Modified 
Gravity may provide a simple, robust, and broadly valid 
benchmark for testing future observations against new 
physics (cf. the role of minimal supergravity for dark 
matter physics). The beyond-Einstein parameterization 
{wo, w a , 7} provides an accessible and reasonable model 
independent approach to studying the physics of the ac- 
celerating universe. 

Acknowledgments 



This work has been supported in part by the Director, 
Office of Science, Department of Energy under grant DE- 
AC02-05CH11231. 



M. Kowalski et al., in preparation 

A. Lue, R. Scoccimarro, G. Starkman, Phys. Rev. D 69, 
124015 (2004) [astro-ph/0401515| 

E.V. Linder, Ph ys. Rev. D 70, 023511 (2004) 
[astro-ph/0402503j 

L. Knox, Y.- S. Song, J.A. Tys on, Phys. Rev. D 74, 
023512 (2006) [astro-ph/0503644| 

M. Ishak, A. Upadhye, D. Sper gel, Phys. Rev. D 74, 
043513 (2006) [astro-ph/0507184] 

P.J.E. Peebles, Large-Scale Structure of the Universe, 
(Princeton U. Press: 1980) 

L. Wang fc P.J. Steinhardt, Ap. J. 508, 483 (1998) 
astro-ph/9804015] 

Rev. 



D 72, 043529 (2005) 



!3.V. Linder, Phys. 

astro- ph/050 7263 

3. Hut erer fc E.V. Lind er, Phys. Rev. D 75, 023519 
(2007) [astro-ph/060868T1 

E.V. Linder, Ph ys. Rev. Lett. 90, 091301 (2003) 
astro-ph/0208512] 

S.V. Linder, Ph ys. Rev. D 73, 063010 (2006) 
astro-ph/0601052] 

C. Wetterich, P hys. Lett. B 594, 17 (2004) 
[astro-ph/0403289j 

M. Doran fc G Robbers, JCAP 0606, 026 (2006) 
[astro-ph/0601544] 

E/Y Linder, Astropart. Phys. 26, 102 (2006) 



a.s tro-ph/0604280 
[15] J.N. Fry, Phys. Lett. B 158, 211 (1985) 
[16] E.V. Linder, Notes on Growth of Structure in a Fried- 

mann Universe, (MPA Internal Research Note: 1988) ; 

E.V. Linder, First Principles of Cosmology, (Addison- 

Wesley: 1997) 
[17] N. Kaiser, MNRAS 227, 1 (1987) 

[18] G. Dvali, G. Gabadadze, M Porrati, Phys. Lett. B 485, 

208 (2000) [hep-th/0005016] 
[19] C. Deffayet, G Dvali, G. Gaba dadze, Phys. Rev. D 65, 

044023 (2002) [astro-ph/010 5068 
[20] F. Perrotta, C. Baccigalupi, S. Matarrese, Phys. Rev. D 

61, 023507 (2000) [astro-ph/9906066] 
[21] C. Baccigalupi et al., in draft 

[22] C. Baccigalupi, S. Matarrese, F. Perrotta, Phys. Rev. D 

62, 123510 (2000) [astro-ph/0005543l 

[23] R. Bean et al Phys. Rev. D 75, 064020 (2007) 

[astro-ph/0611321] 
[24] J.M. Bardeen, Phys. Rev. D 22, 1882 (1980) 
[25] H. Kodama & M. Sasaki, Prog. Th. Phys. Suppl. 78, 1 

(1984) 

[26] CM. Will, Living Rev. Relativity 9, 
[27] L. Pog osian et al., Phys. Rev. D 

[astro-ph/0506396] 
[28] I. Sawicki, Y.- S. Song and W. Hu, Phys. Rev. D 75, 

064002 (2007) astro-ph/0606285 



3 (2006) 

72, 103519 (2005) 



8 



[29] H.F. Stabenau and B. Ja in, Phys. Rev. D 74, 084007 [31] M. Kunz fc D. Sap one, Phys. Rev. Lett. 98, 121301 (2007) 

(2006) [astro-ph/0604038| [astro-ph/0612452] 
[30] E.V. Linder, Jou rnal of Physics A 40, 6697 (2007) 

astro-p h/0610173] 



